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ABSTRACT 

We present observations of GRB 050318 by the Ultra- Violet and Optical 
Telescope (UVOT) on-board the Swift observatory. The data are the first de- 
tections of a Gamma Ray Burst (GRB) afterglow decay by the UVOT instru- 
ment, launched specifically to open a new window on these transient sources. We 
showcase UVOTs ability to provide multi-color photometry and the advantages 
of combining UVOT data with simultaneous and contemporaneous observations 
from the high-energy detectors on the Swift spacecraft. Multiple filters cov- 
ering AA1,800-6,000A reveal a red source with spectral slope steeper than the 
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simultaneous X-ray continuum. Spectral fits indicate that the UVOT colors are 
consistent with dust extinction by systems at z — 1.2037 and z = 1.4436, red- 
shifts where absorption systems have been pre-identified. However, the data can 
be most-easily reproduced with models containing a foreground system of neu- 
tral gas redshifted by z = 2.8 ± 0.3. For both of the above scenarios, spectral 
and decay slopes are, for the most part, consistent with fireball expansion into a 
uniform medium, provided a cooling break occurs between the energy ranges of 
the UVOT and Swift's X-ray instrumentation. 

Subject headings: astrometry - galaxies: distances and redshifts - gamma rays: 
bursts - shock waves - X-rays: individual (GRB 050318) 

1. Introduction 

The multi-instrument Swift observatory (Gehrels et al. 2004) was launched on Nov 20, 
2004. It carries three science instruments, the wide-angle, hard X-ray, Burst Alert Telescope 
(BAT; Barthelmy et al. 2005) which locates GRBs to within 3' on the sky, the narrow-field 
X-Ray Telescope (XRT; Burrows et al. 2005) and the UVOT. Specifications of the UVOT are 
described in Roming et al. (2005). The UVOT instrument has a vital role imaging the field 
containing the burst, minutes after a trigger, and reporting rapidly the afterglow location 
to < 1" accuracy via the GRB Coordinate Network (GCN). UVOTs subsequent role is to 
provide a relatively-uniform sample of the afterglow decay. It is this subsequent role that 
we report on here, describing the first afterglow detected by UVOT in multiple colors, and 
monitoring the decay until 40 ksec after the burst. 

2. Observations 

The Swift-BAT made a 17a detection of GRB 050318 at 15:44:37 UT (Krimm et al. 2005a). 
Burst parameters revised from Krimm et al. (2005c) include a Tgo burst duration of 32 ± 2 
s, with a total fluence of 2.1 xl0~ 6 erg cm -2 in the 15-350 keV band. Within this energy 
band we find evidence of spectral evolution across three peaks in the prompt emission light 
curve. Peak 1, between T-ls and T+5s, (where T is the trigger time), is well fit by a simple 
power law with spectral index (3 BAT = —1.1 ±0.2 (x 2 = 57 for 57 d.o.f.). All uncertainties in 
this paper are reported to a 90% confidence level while spectral and temporal decay indices 
are provided with respect to flux density, e.g., F v ~ t a v" . The burst was quiet for the next 
17 seconds, followed by two overlapping, but resolved peaks. Peak 2 (T+22-27s) fits to a 
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cut-off power law with (3 BAT = -0.2 ± 0.5 and E p = 68±H keV (x 2 = 66 for 56 d.o.f.). The 
spectrum softens considerably during the third peak (T+27-32s), and is fit with (3 BAT = 
0.2 ± 0.4, where E p = 46 ± 7 keV. BAT event data were not recorded for the final 2s of the 
burst. However, examination of the BAT rate data in four energy bands suggests continued 
spectral softening during this period. 

The burst was located to within 3' (90% containment) of RA = 49.651, Dec = —46.392 
(J2000). This corresponds to a Galactic latitude of —55° with a local reddening of E(B- 
V) = 0.018 mag (Schlegel, Finkbeiner & Davis 1998) and a H-equivalent Galactic column 
density of Ah = 2.8 xlO 20 cm" 2 (Dickey & Lockman 1990). After a 54 min delay for Earth 
occultation, Swift slewed so that the narrow-field instruments could monitor the target. 

Within the first 100s settled observation of the UVOT sequence, a V = 17.8 source was 
found 2'. 6 from the BAT position (Fig 1), with no counterpart in archival plates (McGowan 
et al. 2005) and consistent with the ground-based report of Mulchaey & Berger (2005). Sub- 
sequent exposures revealed a fading source, I'M from a transient X-ray counterpart (Mark- 
wardt et al. 2005; Nousek et al. 2005; Beardmore et al. 2005). A complete analysis and 
description of the XRT data reduction is reported in a separate paper (Perri et al. 2005). 
The UVOT position, as reported by De Pasquale et al. (2005), is RA = 03 h 18 m 5 I s . 15, Dec 
= -46° 23' 43".7 (J2000), with 0".3 uncertainties. 

The UVOT completed 36 exposures before GRB 050319 triggered the BAT and be- 
came the new automated target (Krimm et al. 2005b). All detections of GRB 050318 > 2a 
above the background are tabulated in Table 1. AB filter magnitudes and background limits 
are based upon in-orbit zero-point calibrations, and differ from those used by McGowan 
et al. (2005) and De Pasquale et al. (2005), which were based on pre-flight calibrations and 
Vega magnitudes. The afterglow is not detected in UVW1 (centered at approximately 
2,500A), UVM2 (2,200A) or UVW2 (1,800A) light; the first settled exposures yield 3a upper 
limits of 19.3, 19.5 and 21.2 mag, respectively. Detections are made through the U (3,500A), 
B (4.400A) and V (5,300A) filters between T+3,200-5,400s. On the next rotation of the 
filter wheel at ~ T+21,000s the source has decayed below the 3a background threshold in 
both the U and B bands. The magnitudes and detection significances at this epoch are U = 
21.9 ± 0.5, detected at 2.3<r above background and B = 21.6 ± 0.8, 1.3a above background. 
The V source persists for three wheel rotations before fading below the background thresh- 
old between T+23, 000-34, 000s. The U band contains two further marginal source detections 
between 2-3<r above background at T+28,609 and T+40,193. 
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2.1. Source decay 

All U, B and V points, bar the first 2s B exposure, are plotted in Fig. 2. Detections 
> 2a are provided with 90% confidence error bars, while all other points are given as upper 
limits at the 3a level. Using the first three V filter exposures, the powerlaw decay index for 
the V light curve is a v = — 0.87±0.24 (x 2 = 1-1 for 1 dof). Assuming a powerlaw decay, the 
four U band detections, > 2a above background, yield a consistent slope au = —1.00 ± 0.25 
(90% confidence, x 2 — 2 for 2 dof). A weighted- mean of the U and V decay slopes provides 
«u+v = —0.94 ±0.17 and the best fits to U and V data using this slope are plotted in Fig. 2. 
A curve of the same slope is extrapolated to pass through the B detection at T+5,382s. The 
XRT light curve at this epoch, presented by Perri et al. (2005), has a powerlaw decay index 
of 

c^xrt — — 1-2 i 0.1. The probability that gixrt and o;u+v 3xc 

identical is 8%. 

Source detection in the T+21,105s B exposure is significant only to 1.3a and, using the 
B detection at T+5,382s as an anchor point, is inconsistent with the powerlaw decay index of 
a v+v with 99.9% confidence. Either there is under-sampled variability in the source which 
provides us with a biased measure of the decay indices, or we are observing spectral evolution. 
Perhaps it is no coincidence that the next exposure after the second B band observation is 
the one point on the U curve that is an outlier relative to the best-fit powerlaw decay index. 
It is inconsistent with «u+v with 96% confidence. So short-timescale variability is perhaps 
the most plausible interpretation. 



2.2. Spectral properties 

The two absorption systems at z\ = 1.2037 and z 2 = 1.4436 reported by Berger 
et al. (2005) should produce Lyman systems redshifted into the UVM2 band. Assuming 
that consistent non-detections in filters blueward of a particular wavelength reveal either 
dust or the Lyman limit of the host and its redshift, in this section we formally measure the 
spectral slope across the UVOT bands and search for a Lyman edge. All spectral models 
below include Galactic extinction appropriate for the source direction (Sec. 1), using the 
analytic formalism of Pei (1992) with Ry = 3.08. For simplicity, we assume that there is no 
Lya forest in front of the host. 

Using a v+v , Source rates in each filter were interpolated or extrapolated to a common 
epoch of T+4,061s. Using the x 2 fitting method outlined in Arnaud (1996) and references 
therein, a simple powerlaw model to all six points yields a fit with spectral index /3 UV ot — 
—4.9 ±0.5 and x 2 — 24 for 4 dof. The fit is poor statistically and the index is steep compared 
to the 0.2-5 keV slope obtained from simultaneous XRT data of /3 XRT = — l.ll^4 (Sec. 2.3). 
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By adding an absorption edge to the powerlaw model, the quality of the fit improves to a 
statistically acceptable solution; /3 UVOT = —2.4 ±1.5 and, assuming the edge is due to the 
Lyman series, z = 2.8 ±0.2 (x 2 = 1.0 for 3 dof). The best-fit redshift is inconsistent with z 2 
= 1.4436 (from Berger et al. 2005). For comparison, if z is forced to be 1.4436, the best fit 
yields (3 UVO t = —4.5 ±1.5 with x 2 — 15 for 4 dof. For completeness, fixing z at a value of 
1.2037 yields (3 VVOT = -4.5 ± 0.5 (x 2 = 20 for 4 dof). 

We add host extinction to the model above, assuming an SMC grain content with 
Ry = 2.93 (Pei 1992), and coupling the redshift of the dust to the neutral gas. The best 
solution yields /3 UVOT = -1.0±5.0, z = 2.9±0.2 and E(B-V) < 0.26 (x 2 = 1.0 for 2 dof), i.e. 
dust is not a necessary component for this model in order to provide a good description of 
the UVOT data. However, when the dust and gas is re-situated at z 2) the fit does converge 
to an acceptable solution, provided we also include dust and gas at z±, with x 2 — 3.2 for 
2 dof, /3 UV oT = ±1-0 ± 2.0, £(B-V)i = 0.4 ± 0.2 and E(B-V) 2 < 0.27, where £(B-V)i 
and E(B-V)2 are the color excesses at z\ and z 2 respectively. Consequently, gas and dust, 
at the redshifts of the two absorption systems reported by Berger et al. (2005) in front of 
a powerlaw continuum provide an adequate fit to the UVOT data. However, the model 
containing a single gas and dust complex at the larger redshift of z — 2.8 ± 0.2 provides the 
better fit. 

Next we investigate whether the solutions above are biased by assuming an inappropriate 
temporal decay slope during the interpolation of UVOT data to a common epoch. For 
comparison, we repeat the previous exercise using the XRT decay index a XRT = —1.2. Best- 
fit parameters and fit quality vary only a little compared to the previous analysis, and this 
results from the choice of a common epoch which minimizes the systematic uncertainty in 
the interpolation. The best-fit solution without a dust component in the host (which does 
not formally improve the fit) is (3 VVO t = —1-0 ± 5.0 and z = 2.9 ± 0.5 (x 2 = 1.0 for 3 dof). 
The best fit with two dust and gas systems at z 1 and z 2 yields x 2 — 2-6 for 2 dof, j3 VVOT = 
±1.0 ± 4.1, E(B-V)i = 0.5 ± 0.4 and E(B-V) 2 < 0.35. Both solutions are identical to the 
previous analysis within uncertainties. 

Intrinsic continuum slopes in the above models are poorly constrained due to a combi- 
nation of low count rates and the relatively small spectral range of the filters. In the next 
sections, by combining the UVOT data with a simultaneous XRT spectrum, we can place 
further constraints on the UVOT continuum, refine the redshift test and dust measurements, 
and compare a simple fireball model to the data. 
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2.3. Spectral energy density 

Good XRT events have been extracted from within the time interval T+3, 180-5, 822s, 
which is the epoch between the start of the first V exposure and the end of the subsequent 
B detection, and binned by pulse height. The spectral fits below contain a core model of a 
floating powerlaw, combined with fixed quantities for Galactic reddening and extinction in 
the local rest frame (Sec. 1). Galactic abundances are from Anders & Grevesse (1989). We 
use the April 5, 2005 empirical version of the XRT response calibration, which requires an 
additional absorption feature added to spectral models, corresponding to the neutral O K 
feature at 0.54 keV, due to the optical filter. 

The best fit to the XRT data alone yields a powerlaw slope of (3 XKT = —1.2 ±0.3 and an 
integrated 0.2-5 keV flux of (1.7 ± 0.3) xl0~ n erg s" 1 cm -2 . No extra spectral components 
are required with x 2 = 10 for 21 dof. On comparing the UVOT spectrum to the XRT 
spectral model, not only do we find an observed optical/UV spectral index much steeper 
than the X-ray continuum, but also the UV fluxes are > 1 order of magnitude fainter than 
those predicted by the XRT model, therefore the different slopes cannot be caused by a 
spectral break alone. A combined fit of the UVOT and XRT data to the core model yields 
a poor fit with a spectral index (3 = —0.45 ± 0.03 and x 2 — 506 for 28 dof. 

An acceptable combined fit of x 2 — 13 for 25 dof is obtained by adding SMC-like dust 
and neutral gas with Magellanic cloud metallicities (<H/Fe> = —0.5) at one, free floating 
redshift: (3 = -1.0 ±0.1, z = 2.8 ±0.3, E(B-V) = 0.12 ±0.04 and log N H < 2.0xl0 21 cm" 2 . 
The alternative model from Sec. 2.2 replaces the z = 2.8 gas and dust with two systems at 
z 1 and z 2 . Best fit parameters are f3 = -1.1 ± 0.1, E{B-V) 1 = 0.23 ± 0.12, E(B-V) 2 < 0.17, 
N m < 1.7xl0 21 cm" 2 and N m < 1.8 xlO 21 cm" 2 , with x 2 = 16 for 24 dof. These two fits are 
plotted in Fig. 3. Plotting the x 2 landscape of this second model in the E(B-V) 1 -E(B-V) 2 
plane (Fig. 4) reveals that the majority of dust in this scenario is associated with the closer 
of the two systems at z 1 . 

In summary, model fits to the spectrum of GRB050318 at T±4,061s where SMC-like 
neutral gas and dust are situated at z\ and z 2 in front of a simple powerlaw source reproduce 
the UVOT spectral index observed. The spectrum is also well-fit using a single system of 
gas and dust at z — 2.8 ± 0.3. 
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3. Discussion 

3.1. Dust and gas properties 

SMC dust has typically proved to be a good fit to extinction curves in GRB host 
galaxies, e.g. Jakobsson et al. (2004), as one might expect from a host containing a younger 
stellar population (Calzetti et al. 2000). The neutral column density within the host is 
poorly constrained, but by adopting the best fit value from the SMC-dust spectral model 
with z = 2.8, we find a gas-to-dust ratio of N(Hi)/A(V) < 7.6 xlO 21 cm -2 mag" 1 , < 50% 
of typical SMC lines-of-sight (Gordon et al. 2003). If we substitute the SMC dust with the 
Milky Way prescription of Pei (1992; Ry = 3.08) and assume Solar metallicities in the neutral 
gas (Anders & Grevesse 1989) then the fit is also acceptable with x 2 12 for 25 dof, providing 
slightly different best fit parameters of f3 = -1.1 ± 0.1, z = 2.4 ± 0.5, E(B-V) = 0.28 ± 0.13 
and JVh < 1.6 xlO 21 cm" 2 . The gas-to-dust ratio for this case is N(Hl)/A(V) < 3.0 xlO 21 
cm" 2 mag" 1 , consistent with the Galactic mean (Bohlin, Savage & Drake 1978). Similarly, 
LMC dust and gas properties with Ry = 3.16 (also from Pei 1992) yield (3 = —1.0 ± 0.1, 
z = 2.5 ±0.5, E(B-V) = 0.17 ±0.07 and N H < 3.2xl0 21 cm" 2 , X 2 = 13 for 25 dof. In this 
case, N(Hl)/A(V) < 7.5 xlO 21 cm" 2 mag" 1 , which is consistent with the LMC sample from 
Gordon et al. (2003). This general model is therefore acceptable statistically for a range of 
dust and gas content. 

A similar exercise applied to the alternative model, with extinction and absorption 
occurring at Z\ and z 2 , yields a poor fit of x 2 — 50 for 24 dof when Milky Way dust and gas 
populations are assumed. Best fit parameters are (3 = —1.4 ± 0.2, i?(B-V)i = 0.28 ± 0.10, 
E(B-V) 2 = 0.34 ±0.08, N m < 2.0xl0 21 cm" 2 and N U2 < 2.7xl0 21 cm" 2 . An LMC gas 
and dust model results in X 2 = 18 for 24 dof, where (3 = -1.2 ±0.1, E(B-V) 1 = 0.25 ±0.15, 
E(B-V) 2 = 0.17±0.12, JVhi < 3.2xl0 21 cm" 2 and JV H2 < 2.6xl0 21 cm" 2 . The gas-to-dust 
ratio, N(Hi)/A(V), for the LMC model is limited to < 5.0 xlO 21 mag" 1 cm" 2 in the system 
at z 2 , which is a good candidate for the host galaxy; the same ratio for the SMC model is 
unconstrained. While SMC and LMC models provide acceptable fits, the Milky Way model 
does not, although we have made the simplifying assumptions that the dust contents of the 
two systems are identical and that a Lya forest is absent in front of the burst. We also 
note that the constraint from Sec 2.3 that the majority of dust is located in the z = 1.2037 
complex can be dropped if the extinction law in both systems is assumed to be featureless, 
e.g. Savaglio & Fall (2004), where A(X)/R V = E(B - V)(5, 500A/\) S . In this scenario an 
acceptable fit of \ 2 = 16 for 23 dof is obtained where (3 = -1.1 ± 0.1, £(B-V)i < 0.51, 
E(B-V) 2 < 0.46, 5 = 1-6±5:1, N m < 6.4xl0 20 cm" 2 and N m < 4.7xl0 20 cm" 2 . 
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3.2. Interpretation 

The simplest afterglow emission model assumes synchrotron emission from a relativistic 
fireball, expanding into a uniform interstellar medium (Sari, Piran & Narayan 1998; Zhang 
& Meszaros 2004). Assuming that the injection break occurs at an energy < 400 eV then the 
X-ray spectral slope of GRB 050318 indicates that p = 2.4 ±0.2, where (3 = —p/2 according 
to the parameterization of Sari et al. (1998). The simple fireball model then predicts that 
the temporal decay slope should have an index of a = (2 — 3p)/4 = —1.3 ± 0.2. This is 
consistent with the XRT decay index and the B band lower limit, but comparison with the 
other optical bands is less convincing. The U and V indices are consistent with p with 14% 
and 3% confidences respectively. 

If we assume that the cooling break occurs at an energy greater than the injection 
break, and between the UVOT and XRT bandpasses, then the emission models predicts 
c^uvot — 3(1 — p)/A = —1.05 ± 0.2. In this case, the U and V indices are consistent with 
p, with confidences of 80 and 36% respectively. While the best spectral fit from Sec. 2.3 
does not formerly require a cooling break, it does not preclude it either. By replacing the 
model powerlaw continuum from Sec. 2.3 with a broken powerlaw of fixed spectral indices 
Ajvot = —1-05 and /9 XRT = —1.2, an acceptable fit is found using SMC dust with x 2 = 13 for 
25 dof, z = 2.9tol E(B-V) = 0.15 ± 0.02, log N H = 19.9 ± 1.5 cm' 2 and providing a lower 
limit on the cooling break of v c > 4.8 xlO 15 Hz. The model with two dusty sytems at z 1 and 
z 2 yields \ 2 = 16 for 24 dof, £(B-V)i = 0.23 ± 0.12, E(B-V) 2 < 0.17, log N m < 21.1 cm" 2 , 
log iVH2 = 19.7 ± 1.5 cm -2 and v c > 2.4xl0 15 Hz. While the above model of a slow-cooling 
fireball within a uniform ISM fits most of the Swift data well in both redshift scenarios, the 
one caveat is the inconsistency between the B decay index limit and a v+v . 

Assuming z = 2.8 is the host redshift and a cosmological model of Hq = 65 kms -1 
Mpc -1 , Q m = 0.3 and Q\ = 0.7, then the BAT 15-350 keV fluence yields an isotropic energy 
of E iso = (3.6^5'D xlO 52 erg. Transforming a time-averaged value of E p = 49 keV to the rest 
frame of the burst, we find E' p = 1961^ keV (cf E iso = (1.4^; 2 )xl0 52 erg and E' p = 119±|§ 
keV at z 2 = 1.4436). Since both redshifts yield spectral parameters consistent with the 
E iso -E p relationship derived by Amati et al. (2002), the BAT spectral analysis of this burst 
cannot provide useful diagnostics for testing the optical redshift candidates. 

4. Conclusion 

This paper reports the first significant optical detection of a GRB afterglow, and sub- 
sequent monitoring of the decay, by the UVOT instrument on-board the Swift observatory. 
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Compared to a simple powerlaw continuum model, the general deficit of UV emission can 
be fit using either gas extinction at redshifts of z\ = 1.2037 and z 2 = 1.4426, correspond- 
ing to the absorption systems found by Berger et al. (2005), or Lyman depletion from an 
object at z — 2.8 ± 0.3, which would indicate that the two systems at z\ and z 2 belong 
to foreground objects. Consequently the UVOT data cannot unambiguously determine the 
host galaxy redshift for this burst. This will be a common occurence when Swift does not 
detect a UV source. We note that there is no evidence for a z = 2.8 host galaxy in the 
AA5,000-7,000A spectroscopy of Berger et al. (2005), however hosts at this redshift have 
typically been identified by absorption features at bluer wavelengths such as a damped Lja 
line (Hjorth et al. 2003) which would occur at 4,617A at z — 2.8. Since an identified ab- 
sorption line provides only a lower limit to the host redshift, currently available evidence 
only allows a lower limit to be placed on the redshift of z > 1.4436. Decay curves and the 
UVOT/XRT SED reveal mostly consistency with the picture of a slow-cooling fireball in a 
uniform ISM. However the inferred steepness of the B band decay slope, relative to U and 
V, may indicate some deviations from the simple model. 

This work is sponsored at Penn State by NASAs Office of Space Science through contract 
NAS5-00136 and at MSSL and Leicester by funding from PPARC. We gratefully acknowledge 
the contributions of all members of the Swift team. 
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Table 1. UVOT detections (> 2a above background) of GRB 050318 with mid-exposure 

times relative to the trigger (T+), exposure durations, filters, AB magnitudes and 
significance of the detection over background. Filter bandpasses are provided in Roming 

et al. (2005). 
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Fig. 1. — Stacked UVOT-V filter image of the field with the transient source at RA = 03 h 
18 m 51 s . 15, Dec = -46° 23' 43".7 (J2000) and 3' BAT error circle and 6" XRT error circles 
overlaid. Total exposure time for the stacked image is 3,732s. 
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Fig. 2. — U, B and V light curves of GRB 050318. The dashed lines are the best powerlaw fits 
to the U and V time-series, a v+v , excluding upper- limits. The Dotted curve is an identical 
powerlaw model, renormalized to the first-epoch B magnitude. 
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Fig. 3.— The combined UVOT and XRT spectrum of GRB 050318 at epoch T+4,061s, 
compared to two best-fit models, both containing a powerlaw model, reddened and absorbed 
by Galactic material. The Solid line represents this model with an additional system of 
neutral gas and SMC-like dust at z — 2.8. The dotted line represents the model with two 
systems of neutral gas and SMC-like dust at z\ = 1.2037 and Z2 = 1.4436. The dashed line 
is the best-fit intrinsic powerlaw spectrum. 
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Fig. 4. — Confidence map in the J5(B-V)i— J5(B-V)2 plane. The two parameters represent the 
color correction, assuming R v = 2.93, in two SMC-like dusty complexes at z\ = 1.2037 and 
z 2 = 1.4436. Contours are 68, 95 and 99.7% confidence levels and indicate that a significant 
fraction of dust must reside in the closer of the two systems at z\. 



